A comparative study of two phenomenological models of dephasing in series and 

parallel resistors 
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We compare two phenomenological models of dephasing that are in use recently. We show that 
the stochastic absorption model leads to reeisonable dephasing in series (double barrier) and parallel 
(ring) quantum resistors in presence and absence of magnetic flux. For large enough dephasing it 
leads to Ohm's law. On the other hand a random phase based statistical model that uses averaging 
over Gaussian random-phases, picked up by the propagators, leads to several inconsistencies. This 
can be attributed to the failure of this model to dephase interference between complementary electron 
waves each following time-reversed path of the other. 

PACS numbers: 03.65.Yz,73.23.-b,05.60.Gg,11.55.-m 
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I. INTRODUCTION 

Dephasing is defined as the process by which quan- 
tum mechanical interference is destroyed gradually. An 
electron in a sample may lose its phase memory via in- 
teraction with large number of other degrees of freedom, 
like a phonon bath or even due to interaction with all 
other electrons. While the microscopic details of such 
system-bath interactions leading to dephasing is of inter- 
est by itself, we focus on a couple of phenomenological 
models that are currently in use to understand how de- 
phasing affects many quantum interference phenomena 
in mesoscopic systems. In a double slit setup, if ipi and 
'ip2 be the two propagating wave functions that super- 
pose with a phase difference (f>, the interference intensity 
is A = li/^iP + |^/'2p + 2Re[ipiip2 exp(z0)]. In the classical 
limit of complete decoherence the interference term gets 
totally suppressed and ^ = jf/iip -f |'02p- 

Most of the mesoscopic samples have dimensions close 
to the phase coherence length l^, the length scale over 
which the electrons lose their phase memory, can 
be reduced by increasing temperature^. An efficient 
phenomenological model of dephasing was proposed by 
Biittiker— i^. In this model, one attaches 'virtual voltage 
probes' to a system. The probe absorbs phase coherent 
electrons and in turn reinjects incoherent electrons back 
to the system to conserve the overall unitarity of the sys- 
tem plus the probes. This model uses elastic scattering to 
generate dephasing via introduction of the virtual voltage 
probe which carries no net current. A series of side cou- 
pled self-consistent reservoirs, each drawing zero current, 
can effectively induce dephasing; moreover this method 
also allows one to calculate local chemical potentials and 
temperatures*. The main drawback in Biittiker's model 
is that dephasing occurs locally at the point of contact 
between the system and voltage probes. However in a 
natural sample electrons lose phase memory almost uni- 
formly due to interaction with other degrees of freedom. 

This can be taken into account by adding a spatially 



uniform imaginary potential to the Hamiltonian^ one can 
introduce uniform absorption. The main problem with 
this model is that with increase in imaginary potential 
one obtains enhanced back reflection, therefore dephasing 
can not be increased monotonically^'^:^. To avoid this, an 
uniform absorption in the coherent wave function can be 
introduced via a wave attenuation factor—. This factor re- 
duces the wave amplitude by exp(— a^) after traversal of 
a length £ in a free propagating region. Wave attenuation 
added with a proper incoherent reinjection^*' maintaining 
the overall unitarity can be used to introduce dephasing 
in the following way. The three probe Biittiker's model 
can be mapped into an effective two terminal geometry 
by eliminating transmission amplitudes which explicitly 
depend on the third (virtual) voltage probei^. In the 
system, a current I — Ii — ~I2 flows from source (via 
lead 1) to drain (via lead 2). The chemical potential of 
the side coupled voltage probe is adjusted such that no 
current flows through it i.e. ~ 0. The unitarity of 
the 3x3 S'-matrix of the combined system (system plus 
voltage probe) has been used to eliminate the transmis- 
sion coefficients which explicitly depend on the virtual 
voltage probe. Using Landaucr-Biittikcr formul a^'^'^°i^^ 
for coherent transport and eliminating the elements due 
to the third virtual probe, the two probe conductance 
(dimensionless) can be finally expressed asiSii^ 
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The first term is the transmittance T21 from terminal 1 to 
terminal 2. The second term is the incoherent reinjection 
ensuring particle conservation. Rn has the meaning of 
reflectance back into the flrst terminal. All other terms 
in the above expression are self explanatory. All these 
terms are obtainable from a reduced 2x2 non— unitary S- 
matrix which characterizes an effective two terminal ab- 
sorbing device. For further calculation one can use the 



coefficients of an absorbing S'-matrix where absorption 
takes place uniformly via an wave attenuation factor as 
described above. This model of dephasing is known as 
stochastic absorption (SA). For further details one can 
see Ref jiSii^iii^. The SA model is shown to produce rea- 
sonable agreement with experimental results^^. 

In another recently proposed phenomenological 
model'^^ dephasing is introduced via an additive random 
phase (j) of Gaussian distribution to the phase of a 
propagator after it travels a distance £. In this random 
phase (RP) statistical model the following property of 
Gaussian distributed variables is utilized. If the mean 
of the distribution (0) = 0, (exp(i<^)) = exp(— ((/)^)/2). 
Thus if a propagator picks up an extra Gaussian 
random phase (j) of mean zero, it loses an ampli- 
tude exp{—{(fP)/2) in an average sense. This will be 
equivalent to the loss in amplitude in the SA scheme 
[exp(— a£)] if (0^) = 2a£. Utilizing wave interference it 
was 

shown^^ i-jjat; a, Monte-Carlo (MC) averaging over 
many such (perfectly unitary) realizations of random 
phases is expected to generate dephasing. To maintain 
the Onsager-Casimir reciprocity, the unitary S'-matrix of 
the phase randomizer, is symmetric and independent of 
magnetic flujsi^. Thus a propagator picks up the same 
(random) phase in its time-forward and time-reversed 
path. It is possible to do an MC averaging of relevant 
quantities, like conductance, from individual unitary 
processes involving random phase factors to obtain the 
effect of dephasingi^ii^. Other random-phase models, 
like the Lloyd modelii, are also in use. Both the above 
mentioned methods, stochastic absorption and random 
phase model, are expected to allow one to go from 
a fully coherent transport to a fully incoherent one, 
continuously, by varying a from zero to infinity. 

In this paper, we use a series and parallel resistor ge- 
ometry to test these two models for dephasing. We first 
study a quantum double barrier system (Sec. |TT|. The 
two dephasing models predict completely different results 
in the incoherent limit. For the parallel geometry, we 
choose a quantum ring in Scc. lIIII In absence of magnetic 
flux, we study the impact of dephasing on the quantum 
current magnification (CM) ^^i^^ effect. After that, in 
this section, we study the two models of dephasing in the 
context of the Aharonov-Bohm (AB) effect. While RP 
model gives rise to many inconsistencies, SA technique 
produces well behaved predictions. Finally we conclude 
with some discussions in Sec. llVl 



II. DOUBLE BARRIER 

Let us assume a double barrier system consisting of 
two barriers characterized by transmittances Ti and T2. 
We consider the expression for two terminal conductance 
given by G = T or equivalently the two terminal re- 
sistance TZ = 1/T. In the incoherent limit resistances 
should add as in Ohm's law to give the total resistance: 
TZ = 1/G = J2i l/^j ^ system oi i = 1 . . .n barriers 




FIG. 1: Dephasing in quantum double barrier system in 
ID obtained from stochastic absorption model. The data 
is for two rectangular potential barrier having same width 
(wi = W2 = 0.01) and strength {Vi ^ V2 = V = 20). The 
intermediate distance between two barriers is s = 10. We 
plot conductance G (in units of 2e^ /h) as a function of en- 
ergy E/V . The SA gives G = T/2 in the incoherent limit of 
large as, T being the transmittance through a single barrier. 



connected in series. 

To test the SA and RP models of dephasing let us take 
a one dimensional double barrier system consisting of two 
identical rectangular potential barriers of strength V and 
width w. Let the two barriers be separated by a distance 

s. k ^ J2mE/n is the wave- vector and the imaginary 

wave-number is k = j2m{V - E)/fL^. The S'-matrix 
elements for each of the rectangular barriers are. 
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and the S'-matrix is unitary and symmetric. In the SA 
model we assume that the propagator undergoes an at- 
tenuation exp(— as) in a single trip between the barriers. 
Thus the various elements of the combined S-matrix re- 
quired to calculate the conductance in SA method are, 



ti2 — exp(— as 4- iks)/D, 
rii =r + r cxp(-2as + 2iks) /D 



(3) 



In the above expressions, D 
Using these and Eq.lJT]) we find 



with t2l tl2, 7-22 = rii 

1 — exp(— 2q;s -I- 2iks). 
a loss of oscillation (coherence) in the total conductance 
along with a decay in total conductance as a function 
of a (See Fig[T|). At large enough absorption as ^ 00 
(practically already for as > 5) the total conductance 
asymptotically approaches the value G T/2, a result 
expected from the earlier discussions on incoherent limit 
{1/G = 1/T + 1/T). 

For the implementation of MC averaging over random 
phases, note that the transmittance through a single such 
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FIG. 2: Dephasing in quantum double barrier system in ID 
using random phase model. The system parameters are the 
same as in Fig[T] We plot conductance G (in units of 2e^ /h) 
as a function of energy E /V . The MC simulation predicts 
G = T/{2 — T) in the incoherent limit of large as, T being 
the transmittance through a single barrier. 



barrier is T = + (fc^ + k^)^ smh^ {2kw) / 'ik'^ k^] and 
reflectance \s R = \ — T. We assume that the propaga- 
tor picks up a random phase while traversing the free 
space (s) between the barriers, the total transmittance 
through the double barrier system. For each realization 
of the random phase, total transmittance across the dou- 
ble barrier system is 



'12 



l-2Rcos{2ks + 2(j))+R'^' 



(4) 



The random phase is assumed to follow a Gaussian dis- 
tribution of mean (0) = and variance {(jp) = 2as. An 
average over 500 realizations of the random phase is per- 
formed at each a to observe a monotonic decay of the 
oscillations (dephasing) in total conductance G = {Ti2). 
In the limit of large dephasing factor as — s- oo (prac- 
tically aheady for as > 1) G ^ T/{2 - T) (FigH]). 
This is what one obtains by adding the S'-matrices of 
the two barriers incoherentljsiiii^. For two dissimilar 
barriers the RP method in the incoherent limit predicts 
G = TiT2/(l-i?ii?2), G = r/(2-T) being a special case 
with Ti = T2 = T. This result is in clear disagreement 
with the expectation of Ohm's law in series resistors in 
the limit of complete dephasing. 



III. QUANTUM RING 

We consider a quantum ring (Figl3]) where L„, Ld are 
the lengths of the upper and lower arms respectively. The 
total circumference of the ring is L = L„ + Ld- The ring 
is connected to two leads at junctions Ji and J2. These 
connections are described by a scattering matrix Sc- For 
the junction Ji the outgoing amplitudes {a',b'i,c'i) are 
connected to the incoming amplitudes (a, 6i,ci) via Sc- 




FIG. 3: A schematic diagram of the ID quantum ring consid- 
ered. $ denotes the flux enclosed. The arrows show counter 
propagating amplitudes of wave functions. fii^2 chemical po- 
tentials at two reservoirs connecting the ring at the junctions 
Ji,2 by two leads. 
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FIG. 4: Results from stochastic absorption in quantum ring 
with Lu/L = 0.45 and Ld/L = 0.55 in absence of magnetic 
flux. The coupler e — 4/9. Conductance G (in units of 2e^ /h) 
and circulating current density jc (in units of ji„, see text) 
in a open quantum ring with changing Fermi momentum kL. 
The amount of circulating current jc and the range of energy 
over which it is obtained reduces with increase in dephasing 
a. At large enough a{— 5) G ^ e. 



Similarly for junction J2, {d' h'^ = Sc{d 62 . 
junction S-matrices can be modeled as^S 



Sc 




The 



(5) 



with p = y\ - 2e- l]/2 and q = [^1 - 2e+ 1]/2. The S- 
matrix Sc is characterized by a single parameter e, which 
can take any value from (decoupled) to 0.5 (strongly 
coupled). The quantum mechanical continuity of wave- 
function and its derivative is attained at the junctions 
for e = 4/9. Without any lack of generality we use this 
value for the coupler e in our study. 

The amplitudes 61 , 62 and b'^ , b'2 of the upper arm are 
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FIG. 5: Results from random phase model in the quantum 
ring. Conductance G (in units of 2e'^ /h) and circulating cur- 
rent density jc (in units of jin) in a open quantum ring as a 
function of Fermi momentum kL. With increase in a amount 
of Jc decreases. Notice the occurrence of jc 7^ for a — 0.05 at 
some Fermi momenta (e.g. as denoted by the arrow) although 
for Q = at the same momenta jc = 0. At the incoherent limit 
of large a{= 5) conductance obtained from MC averaging is 
G < e. The parameters used are Lu/L = 0.45, Ld/L — 0.55 
and e = 4/9. 



related by an S'-matrix 



e.(fe-e/L)L„p^ 
^i(k+e / L) Q 



Similarly for lower arm 



^i{k-0/L)Ld p Q 



.(6) 



■(7) 



Here kLu {k Ld) denotes the phase picked up by the 
electron while traversing the upper (lower) arm of the 
ring in absence of magnetic flux and any dephasing fac- 
tor. In presence of Aharonov-Bohm flux ($) phase of 
the wave-function in upper (lower) arm gets shifted by 
an amount OL^/L (OLd/L). This phase-shift in two 
arms add up to give the total flux piercing the ring i.e. 
9 Lu/L + 9 Ld/L = 27r$/$o, where $0 = hc/e, flux 
quanta. An extra factor P„ {Pa) in the propagator along 
the upper (lower) arm introduces dephasing in the sys- 
tem. 

Before going into the calculation of conductance and 
impact of dephasing due to the two dephasing-models we 
study, let us first discuss what we expect in the com- 
pletely incoherent limit. The incoming beam transits 
electrons in the upper (lower) arm with transmission co- 
efficient e. Hence the upper (lower) arm have same resis- 
tances at the two junctions 1/e. In the incoherent limit 
the resistances in upper (lower) arm add as in Ohm's law, 
TZu = 2/e. Now these resistances add in parallel leading 
to the total conductance G = 2/TZu = e. 

In the method of stochastic absorption Pu = 
exp(— aLtj) and Pd = exp(— aLd) act as the continu- 



FIG. 6: Stochastic absorption results in Aharonov-Bohm os- 
cillations. The ring geometry is L^/L — 0.45, Ld/L = 0.55. 
The coupler strength e = 4/9. We fixed kL = 5. The visibil- 
ity of AB oscillations in G (in units of 2e^ /h) clearly dies out 
with increasing a to asymptotically reach a fiux independent 
conductance G = e in the limit of large a(= 5). 



ous lossy channels leading to dephasing. Thus combin- 
ing these iS-matrices one can obtain the elements of the 
effective two terminal 'non-unitary' S'-matrix and using 
Eq.([T]) find out the conductance G in presence of lossy 
channels. In absence of magnetic flux {9 = 0) we study 
the conductance across an asymmetric ring {Lu 7^ Ld). 
Throughout this analysis we use L„ = 0.45 and Ld — 0.55 
in units oi L = L^ + Ld = 1. The conductance in such 
an asymmetric ring is known to show Fano-resonances as 
well as Breit-Wigner line shapes as a function of Fermi 
energjiiSii^i^. The Fano resonances accompany occur- 
rences of circulating current in the ring. In presence of 
transport current through an asymmetric ring system, 
depending on Fermi-energy, in one of the arms current 
can be larger than the transport current. Thus it is 
known as current magnificat io n^^i^^'^^ . In such a case, 
to maintain current conservation, in the other arm, cur- 
rent flows opposite to the bias. The amount of current 
flowing opposite to the bias (negative current) in upper 
and lower arms of the ring gives the magnitude of circu- 
lating current density j^. We assign positive (negative) 
sign to jc for clockwise (anti-clockwise) circulating cur- 
rent density i.e. negative current density in lower (upper) 
arm jd < {ju < 0). We measure the current densities 
in the unit of incident current density jin. In a small 
energy interval dE about the Fermi energy, the incident 
current density is jin — ev-^f{E)dE, where f{E) is the 
Fermi distribution function, ^ = 2/hv is the density of 
states (DOS) in the perfect wire and v = Tik/m. For 
the zero temperature calculations f{E) — 1 for occupied 
states. Thus the incident current density becomes jin — 
(2e/h)dE. In dimensionless units the current density in 



upper (lower) arm is ju/jin = \b'i\^ - 

(Wi»=Kp-|Cip = |c2| 




FIG. 7: Results of random phase model in Aharonov-Bohm 
oscillations. The parameters used are Lu/L = 0.45, Ld/L = 
0.55, kL — 5 and e — 4/9. The method clearly fails to reduce 
the visibility of AB oscillations in G (in units of 2e^ /h) with 
increasing a beyond 5. 



Occurrence of circulating current is a purely quantum 
phase coherence effect. Thus with increasing dephasing 
the circulating current density is expected to decay, van- 
ishing completely for fully incoherent transport. So we 
focus on the change in this quantity with increasing a to 
quantify the degree of dephasing introduced by the two 
dephasing models at hand. FigH] shows the total con- 
ductance G and the circulating current density jc as a 
function of Fermi momentum {kL) in absence (a = 0) 
and presence {a = 0.1) of dephasing. It is clear from the 
figure that amount of the circulating current density jc 
and the range of Fermi momentum over which it is found, 
decreases with increasing a. It is assumed here that the 
reinjected current which leads to classical behaviour does 
not contribute to circulating current. At the a — > oo limit 
it is found that G ^ e [see C'ia = 5) in FigH]. In this 
limit the circulating current density is completely absent. 
Thus predictions from SA is completely consistent with 
our expectations about the incoherent limit. 

In the RP model, on the other hand, a random phase 
is picked up by the propagator while traversing the two 
arms of the ring. Dephasing in this model is introduced 
by Pu = exp(i0„) and Pd = exp(i(/)d) such that the 
phases obey Gaussian distributions with zero mean and 
variance cr^ = 2aLu and = 2aLd- At each realiza- 
tion of the random phases the S-matrices connecting the 
left and right junctions (Eqs. (O and ([7])) are unitary 
and symmetric. Thus they obey Onsager reciprocity^^ 
at each given realization. In this method one needs to 
find the transmittance T12 in each realization and per- 
form an MC averaging over various realizations. We av- 
eraged over 500 realizations for each value of a. From 
Figl5]it is clear that this method is also capable of reduc- 
ing the amount of circulating current density for most of 
the Fermi momenta. However, at some regimes of Fermi 
momenta, e.g. near kL = 5.6, we now observe non-zero 
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FIG. 8: First three harmonics at, i — 1, 2, 3 of G as a function 
of dephasing factor a. Result from stochastic absorption. 

jc in presence of non-zero q;(=0.05) though at the same 
points with a = circulating current density was ab- 
sent (Figini). This clearly goes against the basic notion 
of dephasing ~ as this is enhancement of an effect (cir- 
culating current) purely attributed to phase coherence^. 
SA method is seen to be free from such discrepancies. 
However at large enough a(=5) the circulating current 
density completely vanishes and the total conductance 
becomes independent of kL. The a — s- cx) limit of G ob- 
tained by this method remains slightly smaller than that 
obtained in SA method (G < e, see Figl5]). 

Next we switch on the magnetic fiux $ and study how 
dephasing affects the Aharonov-Bohm (AB) oscillations 
of conductance G as a function of magnetic flux. The 
SA absorption method leads to a monotonic decay in the 
visibility of oscillations (difference between maxima and 
minima) with increasing a (Fig[6|) eventually the conduc- 
tance becoming independent of $ in the limit aL — s- 00. 
In fact at aL = 5 it reaches the incoherent limit of G — > e. 
On the other hand, though the RP model leads to de- 
phasing at small a, it rapidly saturates and fails to kill 
the AB oscillation with increasing a any further. 

To understand the results, we explicitly extract the 
various harmonics in the conductance data presented 
above. The n-th harmonic is 

a„ = - / Gcos{n9)d9. (8) 
TT Jo 

The SA method shows fast decay in the various harmon- 
ics. In FiglHlit is clearly seen that the higher is the har- 
monic the faster it decays to zero with increasing a. How- 
ever, the same first three harmonics from MC data shows 
a different behavior (Figl5]). All the odd harmonics de- 
cay to zero, higher harmonics going to zero faster than 
the lower one. Whereas the even harmonics show an ini- 
tial decay followed by a saturation. It is clear from Fig[9] 
that the second harmonic a2 saturates to 02 = —0.14 
with increasing a. Therefore the AB oscillation that was 
predominantly $0 oscillation at aL = becomes pre- 
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FIG. 9: First three harmonics at, i — 1, 2, 3 of G as a function 
of dephasing factor a. Result from random phase model. 

dominantly a $0/2 oscillation for aL > 5 (Figl?]). In 
fact, the RP model fails to kill all the even harmonics in 
the AB oscillations. 

This failure of RP model can be understood in the fol- 
lowing way. Consider two paths which are time reversed 
version of each other. Start from a point to go around 
the ring (clockwise) once and come back to origin. One 
can also travel anticlockwise and come back to the same 
point. These are the two time reversed paths. These two 
paths together encloses the flux twice. This leads to the 
$0/2 oscillation as a function of magnetic flux. Notice 
that the phase difference, due to the magnetic flux, be- 
tween these two paths is 47r$/$o- However, due to the 
symmetric nature of the phase randomizing 5'-matrix, 
each of these two paths picks up the same random phase. 
Thus in the phase difference, the contribution from ran- 
dom phase cancels out. Therefore phase randomization 
in the RP model fails to kill even periodicity ($o/2n with 
integer n) contributions in one dimensional geometries. 
This is true for all even harmonics. 



IV. DISCUSSIONS AND CONCLUSION 

Studies on series and parallel geometry showed that 
the SA method leads to monotonic dephasing (decay in 
visibility of oscillations or purely coherent effect like cir- 
culating current densities in the ring) with increase in 
dephasing parameter a. In the fully incoherent limit of 
a ^ 00 we recover Ohm's law, inverse transmittance 
through individual elements (a barrier in double barrier 
case and an arm in quantum ring) behaving as classical 
resistances. On the other hand RP model of introduc- 
ing random phases to the propagator and averaging over 
many realizations leads to several inconsistencies. For 
double barrier transmittance, this leads to an incoherent 
limit that differs from the classical Ohm's law. For the 
transmission through quantum ring, this method gener- 
ates circulating current in some regimes of Fermi wave- 
vectors where there was no circulating current in absence 
of phase randomization. This behavior is clearly beyond 
the very principle of dephasing, as the circulating cur- 
rent is of purely quantum mechanical origin. If the ring 
is penetrated by magnetic flux, RP model fails to kill 
the $0/2 oscillations. This is because RP model fails 
to dephase interference between complementary electron 
waves which follow time reversed path of each other. 

Thus in conclusion, we find the SA method to be a 
reliable phenomenological technique of introducing de- 
phasing. This systematically increases dephasing with 
increasing value of the dephasing parameter a, finally 
leading to classical Ohm's law in the limit of large a. 
However the RP model fails in many aspects to be a use- 
ful phenomenological model of dephasing. 
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